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Momentum correlation functions of light (anti)nuclei formed by the coalescence mechanism of 
(anti)nucleons are calculated in several relativistic heavy-ion collision systems, namely i?B +5° B, 
160 +2° O, 38Ca +59 Ca as well as 73’Au +72" Au, at different centralities at center of mass energy 
V5NN = 39 GeV within the framework of A Multi-Phase Transport (AMPT) model complemented 
by the Lednicky and Lyuboshitz analytical method. Momentum correlation functions for identical 
or nonidentical light (anti)nuclei are constructed in the above collision systems at so high collision 
energy. The results suggest that emission of light (anti)nuclei occurs from a source of smaller space- 
time extent in more peripheral collisions. The effect of system-size on the momentum correlation 
functions of identical or nonidentical light (anti)nuclei is also explored in several central collisions. 
The results indicate that the emission source-size of light (anti)nuclei pairs deduced from their 
momentum correlation functions and system-size is self-consistent. Momentum correlation functions 
of nonidentical light nuclei pairs gated on velocity are applied to infer the average emission sequence 
of them. The results indicate that protons are emitted in average on a similar time scale with 
neutrons but earlier than deuterons or tritons in the small relative momentum region. In addition, 
larger interval of the average emission order among them is observed by large centrality and smaller 


system collisions. 


I INTRODUCTION 


In heavy-ion collisions (HICs), two-particle momentum 
correlation function is different from the original applica- 
tion in astronomy [1, 2], and has been normally utilized 
to extract space-time information of the emission source 
and probe the dynamical evolution of nuclear collisions 
in an extensive energy range [3-12]. Many different stud- 
ies on the two-particle momentum correlation functions 
in intermediate energy HICs can be also found in liter- 
ature, eg. Refs. |11, 13-25], which include the momen- 
tum correlation functions of neutron, proton as well as 
light charged particle (LCP) pairs. Multi-variable de- 
pendences of the momentum correlation functions, such 
as impact parameters, total momentum of particle pairs, 
isospin of the emission source, nuclear symmetry energy, 
nuclear equation of state (EOS) as well as in-medium 
nucleon-nucleon cross section (NNCS) etc., contain a 
wealth of information about the space-time character- 
istics of intermediate energy HICs. In high energy HICs, 
two-hadron momentum correlation functions, also called 
as Hanbury Brown-Twiss (HBT) interferometry, were 
also well measured and some interesting properties on 
emission source were extracted [26, 27]. Oscillations of 
the extracted HBT radii versus emission angle indicate 
sources elongated perpendicular to the reaction plane. 
The results indicate that initial shape could be remained 
to be identified even though the pressure and expansion 
the collision system. Furthermore, interaction between 
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antiprotons has been also measured with the momen- 
tum correlation functions and the equality of interac- 
tions between p-p and p-p was approved by the STAR 
Collaboration [28]. The interaction property of the par- 
ticle pairs has been discussed for other particles, for in- 
stance A pairs [29], proton-Q and proton-= etc [30, 31], 
with the same momentum correlation technique. Fur- 
thermore, the measurements of momentum correlation 
functions for nonidentical nucleons and light clusters can 
be used to characterize the mean emission sequence of 
them, which was firstly proposed in Ref. [32]. Theoretical 
study has been extended to different kinds of nonidenti- 
cal particle pairs, for instance p-d, n-p [33-36], «-p [37], 
K+-K- [38], d-t [11, 22] as well as ?He-a particles [39] 
in intermediate energy HICs. 

In this work we extend, for the first time, the momen- 
tum correlation functions of light (anti)nuclei to ultra- 
relativistic heavy-ion collisions simulated by A Multi- 
Phase Transport (AMPT) model [40, 41] coupled with 
a dynamical coalescence model [42-44], specifically at 
/snn = 39 GeV. Different gating conditions such as cen- 
trality gates, system-size gates as well as velocity gates 
are applied to the momentum correlation functions of 
light (anti)nuclei pairs. In particular, we report on the in- 
dication of the emission chronology of protons, deuterons 
and tritons which can be deduced from their correspond- 
ing momentum correlation functions in ultra-relativistic 
HICs at J/sww = 39 GeV. The emission sequence inferred 
from the correlation functions is expected measurable in 
future experiments to verify our deduction from the co- 
alescence picture. 


The rest of this article is organized as follows. In Sec- 
tion II A and II B, we briefly describe A Multi-Phase 


Transport model |40, 41] and the coalescence model |42— 
44], then introduce how to calculate the momentum cor- 
relation functions of particle pairs by using the Lednicky 
and Lyuboshitz analytical formalism [3, 45—48| in Sec- 
tion II C. In Section III, we summarize the simulated re- 
sults of the light (anti)nuclei momentum correlation func- 
tions gated on various parameters in relativistic heavy- 
ion collisions. Section III A compares the results of 
proton-proton and proton-antiproton momentum corre- 
lation functions with experimental data from the RHIC- 
STAR collaboration. From Section III B to III D, iden- 
tical and nonidentical light (anti)nuclei momentum cor- 
relation functions gated on different conditions are sys- 
tematically discussed, respectively. Finally, a summary 
and outlook are given in Section IV. 


II. MODELS AND FORMALISM 
A. AMPT model 


To obtain phase-space distributions of (anti)particles, 
A Multi-Phase Transport model [40, 41] is used as the 
event generator, which has been applied successfully for 
studying heavy-ion collisions at relativistic energies, eg. 
[43, 44, 49-57]. We briefly review the main components 
of the AMPT model used in the present work. In the 
version of AMPT, the initial phase-space information of 
partons is generated by the heavy-ion jet interaction gen- 
erator (HIJING) model [58, 59]. The interaction between 
partons is then simulated by Zhang's parton cascade 
(ZPC) model [60]. During the hadronization process, 
a quark coalescence model is used to combine partons 
into hadrons [61-63]. Then, the hadronic rescattering 
evolution is described by a relativistic transport (ART) 
model [64]. 


In this paper, the collisions of 19B +£° B, 190 +46 O, 
30Ca +59 Ca for the 0 — 10 % central collisions at mid- 
rapidity (|y| < 0.5) as well as 137 Au +42” Au at same 
mid-rapidity for five centralities of 0 — 10 96, 10 — 20 
70, 20 — 40 %, 40 — 60 %, and 60 — 80 % at SNN 
— 39 GeV are simulated. The phase-space distribu- 
tions of (anti)particles are selected at the final stage in 
the hadronic rescattering process (ART model [64]) with 
considering baryon-baryon, baryon-meson, and meson- 
meson elastic and inelastic scatterings, as well as reso- 
nance decay or week decay. The transverse momentum 
spectra of light (anti)nuclei have been successfully repro- 
duced by the AMPT model with the maximum hadronic 
rescattering time of 100 fm/c [44]. Therefore, the same 
maximum hadronic rescattering time is used for the most 
calculations in this work except for a quantitative com- 
parison with a preliminary p-p and p-p data from the 
STAR collaboration in Sec. III A. 


B. Coalescence model 


In our model calculations, light (anti)clusters such as 
(anti)deuterons and tritons are constructed by using the 
coalescence model as follows [65, 66]. The coalescence 
model has been used widely in describing the production 
of light clusters in the intermediate [67-71] and high ener- 
gies [72, 73] heavy-ion collisions. In the model, the prob- 
ability for producing M-nucleon cluster is determined by 
its Wigner phase-space density and the nucleon phase- 
space distribution at the freeze-out stage [42]. The mul- 
tiplicity of an M-nucleon cluster in transport model sim- 
ulations for heavy-ion collisions is given by, 
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where T;,,T;, and ki, kiu. , are the relative coordi- 
nates and momentum in the M-nucleon rest frame, and 
spin-isospin statistical factor G is 3/8 for (anti)deuteron 
and 1/3 for triton [42]. In addition, p" is the Wigner 
density function, which is different for all kinds of parti- 
cles. Therefore, we will calculate separately the Wigner 
phase-space density of (anti)deuteron and triton in de- 
tail. The Wigner phase-space density of (anti)deuteron 
is constructed by, 
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where k — (kı 
r = (ri — T3) is the relative coordinate of (anti)proton 
and (anti)neutron. The Wigner phase-space density of 
triton is constructed by a spherical harmonic oscilla- 


tor [42, 43, 74], 


— ka) /2 is the relative momentum and 
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where p and A are relative coordinates, ko and ky are the 
relative momenta in the Jacobi coordinate. 

The above parameters of the Gaussian fit coefficient 
c; and w; for (anti)deuteron as well as b for triton are 
given in Ref. [42]. Based on the phase-space informa- 
tion of light (anti)cluster obtained by the above coa- 
lescence model, the momentum correlation functions of 
(non)identical light (anti)cluster pairs can be discussed 
in the following of this paper. 


C. Lednikynd Lyuboshitz technique 


Finally, we briefly review the technique of the two- 
particle momentum correlation function proposed by 
Lednický and Lyuboshitz [45-47]. The method calculat- 
ing the two-particle momentum correlation function in 
heavy-ion collisions is based on the principle as follows: 
when two particles are emitted at small relative momen- 
tum, their momentum correlation function is determined 
by the space-time characteristics of the production pro- 
cesses owing to the effects of quantum statistics and final- 
state interactions [3, 48]. Therefore, the two-particle mo- 
mentum correlation function can be expressed through a 
square of the symmetrizied Bethe-Salpeter amplitude av- 
eraging over the four coordinates of the emitted particles 
and the total spin of the two-particle system, which rep- 
resents the continuous spectrum of the two-particle state. 
In the method, the final-state interaction of the particle 
pairs is assumed independent in the production process. 
According to the conditions in Ref. [33], the momentum 
correlation function of two particles can be written as the 
expression: 
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where r* = x4 — X2 is the relative distance of the two 
particles in the pair rest frame at their kinetic freeze- 
out, k“ is half of the relative momentum between two 
particles in the pair rest frame, S (r*, k^) is the proba- 
bility to emit a particle pair with given r* and Kk, i.e., 
the source emission function, and Yy» (r*) is the reduced 
Bethe-Salpeter amplitude which can be approximated by 
the outer solution of the scattering problem |28, 75]. This 
approximation is valid on condition |t*| < m (r*)^, which 
is well fulfilled for sufficiently heavy particles like protons 
or kaons and reasonably fulfilled even for pions [46]. In 
the above limit, the asymptotic solution of the wave func- 
tion of the two charged particles approximately takes the 
expression: 


Wy (r*) = e?«4/ A, (A) x 


e p CA LE) + für) £M . (5) 


r* 


In the above equation, 6. = arg I (1 + iA) is the Coulomb 
s-wave phase shift with \ = (k*a;) ! where a, is the two- 
particle Bohr radius, Ae (A) = 2zxA [exp (2A) — 1| is 
the Coulomb penetration factor, and its positive (neg- 
ative) value corresponds to the repulsion (attraction). 


G (p, À) = 4/ Ac (A) [Go (p, A) + iFo (p, A)| is a combina- 
tion of regular (Fo) and singular (Go) s-wave Coulomb 
functions [46, 47]. F (—ià,1,i£) = 1 + (—iA) (i£) /1!? + 
(—iX) (—iA + 1) (i£)? /2!? +--+ is the confluent hyperge- 
ometric function with € = k*r* + p, p = k*r*. 
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fc (k*) is the s-wave scattering amplitude renormalizied 
by the long-range Coulomb interaction, with h(A) = 
A? Sor, [n (n? + X2)] ^ __C—In [A] where C = 0.5772 is 
the Euler constant. Ke (k*) = m + idok” PEP dave is 
the effective range function, where dọ is the effective ra- 
dius of the strong interaction, fo is the scattering length 
and P is the shape parameter. The parameters of the 
effective range function are important parameters char- 
acterizing the essential properties of the final-state inter- 
actions, and can be extracted from the correlation func- 
tion measured experimentally [28, 34, 76, 77]. Table I 
shows the parameters of the effective range function for 
different particle pairs in the present work. 


TABLE I. Experimental determination of the effective range 
function parameters for n-n (n-n), p-p (p-p), t-t, n-p (n-p), 
p-d (p-d), p-t and d-t systems [28, 76, 77]. 


System Spin fo (fm) do (fm) P (fm? 


n-n (n-n) 0 17 2.7 0.0 
p-p (p-p) 0 7.8 2.77 0.0 
it 0 1x10? 0.0 0.0 
n-p (n-p) 0 23.7 2T 0.0 
p-d (pa) 1/2 -2.73 207 0.08 
3/2 -1L88 2.63 -0.54 

p-t 0 1x10? 0.0 0.0 
d-t 0 1x10? 0.0 0.0 


In the above table, for n-n (n-n) and n-p (n-p) mo- 
mentum correlation functions which include uncharged 
particle, the Coulomb penetration factor (A, (A)) is not 
considered and only the short-range particle interac- 
tion works. For the momentum correlation functions 
of charged particles such as p-p, p-p (p-p), d-d (d-d), 
t-t, p-d (p-d), p-t and d-t, both the Coulomb interac- 
tion and the short-range interaction dominated by the 
s-wave interaction are taken into account. The momen- 
tum correlation function of p-p (p-p) particle pairs is 


dominately contributed by only the singlet (S — 0) s- 
wave final-state interactions while both spins 1/2 and 
3/2 contribute in the case of p-d (p-d) system. Moreover, 
for (anti)deuteron-(anti)deuteron momentum correlation 
function, a parametrization of the s-wave phase shifts 6 
has been used from the solution of K, (k*) = cotó for 
each total pair spin S = 0,1,2. Note that the effective 
range function for the total spin S = 1 is irrelevant, since 
it does not contribute due to the quantum statistics sym- 
metrization. 


III. ANALYSIS AND DISCUSSION 


A. Comparison between our simulations and 
preliminary data 


Fig. 1 presents results of p-p and p-p correlation func- 
tions for three different centrality classes of 0 — 10 96, 
10 — 30 96, and 30 — 70 96 calculated by the AMPT model 
in Au + Au collisions at J/syw = 39 GeV. Within the 
cut of transverse momentum p; and rapidity y, we con- 
front the experimental data with the predictions of the 
AMPT model combined with Lednicky and Lyuboshitz 
code. When the phase-space information of nucleons at 
the maximum hadronic rescattering time of 700 fm/c 
is selected from the AMPT model, it is found that the 
results can well describe the experimental data for the 
momentum correlation functions of p-p and p-p from the 
RHIC-STAR collaboration [78, 79], especially in more 
central collisions. Considering that the preliminary ex- 
perimental results were not corrected by feed-down ef- 
fect corrections |78, 79], the real correlation functions 
for primary p-p and p-p could be much more stronger. 
In this case, using much longer hadronic rescattering 
time of 700 fm/c in the AMPT model might be a rea- 
sonable choice for making quantitative comparison with 
feed-down uncorrected data since the system will become 
more expanded and weakly correlated among particles 
after longer hadronic rescattering time in AMPT. How- 
ever, the quantitative reproduction is not our main con- 
cern in the present work. In the following calculations, 
we fixed the maximum hadronic rescattering time at 100 
fm/c and presented systematic results among different 
light (anti)nuclei. 


B. Centrality and system-size dependence of 
identical light (anti)nuclei momentum correlation 
functions 


The centrality dependence of the two-particle mo- 
mentum correlation function can systematically investi- 
gate the contributions from the system-size and parti- 
cle interactions on the correlations. Fig. 2 (a) and (c) 
presents the momentum correlation functions of iden- 
tical (anti)particle pairs (n-n (n-n) and p-p (p-p) ) for 
197 Au +29" Au collisions at different centralities of 0 — 10 
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FIG. 1. Proton-proton (a) and proton-antiproton (b) mo- 


mentum correlation functions for different centrality classes in 
v5NN = 39 GeV Au + Au collisions. Solid markers represent 
the preliminary experimental data from the RHIC-STAR col- 
laboration [78, 79], and lines represent our model calculation 
results from the AMPT model plus the Lednicky and Lyu- 
boshitz code. Note that the longer hadronic rescattering time 
of 700 fm/c is used in this specific calculation in comparison 
with the data. 


76, 10 — 20 96, 20 — 40 96, 40 — 60 96, and 60 — 80 % at 
V'5NN = 39 GeV. The momentum correlation functions 
of (anti)neutron pairs exhibit more than unity in Fig. 2 
(a), which is caused by the attractive s-wave interaction 
between the two (anti)neutrons. In Fig. 2 (c), the shape 
of the (anti)proton—(anti)proton momentum correlation 
functions looks as expected from the interplay between 
the quantum statistical (QS) and final state interactions 
(FSI) and is consistent with previous results [13, 28, 34]. 
The (anti)proton—(anti)proton momentum correlation 
functions exhibit less than unity at low relative mo- 
mentum q in Fig. 2 (c), which is mainly caused by the 
Coulomb repulsion between the (anti)proton pairs. With 
increasing relative momentum, the attractive s-wave in- 
teraction between the two (anti)protons gives rise to a 
maximum of the (anti)proton—(anti)proton momentum 
correlation functions at q around 0.020 GeV in Fig. 2 
(c). The antiproton—antiproton momentum correlation 
functions show a similar structure with proton pairs, re- 
sulting from the same attractive interaction between two 
antiprotons [28]. Fig. 2 (a) and (c) compare five central- 
ities of 0 — 10 96, 10 — 20 96, 20 — 40 96, 40 — 60 %, and 
60 — 80 % of the two-(anti)particle momentum correla- 
tion functions. The enhancement strength of the n-n (n- 
n) and p-p (p-p) momentum correlation functions is ob- 
served in peripheral selection. These results indicate that 
(anti)particle emission occurs from a source of smaller 
space—time extent in peripheral collision. In addition, 
the effect of system—size on the momentum correlation 
functions of (anti)particles is also investigated by four dif- 
ferent systems, namely 19 B +Ł0 B, 10 4-16 O, 38Ca +38 Ca 


and 19" Au 4-2" Au, in central collisions. In Fig. 2 (b) and 
(d), the n-n (n-n) and p-p (p-p) momentum correlation 
functions appear very sensitive to system-size and an en- 
hancement strength is observed when particle pairs emit- 
ted from smaller system collisions. This enhancement 
strength of the momentum correlation functions for par- 
ticle pairs is a physical effect stemming from the smaller 
space—time extent of the emission source [8]. Therefore, 
the emission source-size of particle pairs obtained by their 
momentum correlation functions and system-size is self- 
consistent. 
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FIG. 2. The momentum correlation functions at mid- 


rapidity (|y| « 0.5) of (anti)neutron-pairs and (anti)proton- 
pairs as a function of five different centralities for 19" Au 4-19" 
Au reaction at J/sww = 39 GeV are presented in (a) and 
(c), respectively. The momentum correlation functions of 
(anti)neutron-pairs and (anti)proton-pairs at mid-rapidity 
(|y| < 0.5) for 0-10% central collisions of §°B+3°B, 10--1* O, 
20Ca +39 Ca as well as 19" Au 4-12" Au systems at V SNN = 39 
GeV are presented in (b) and (d), respectively. The p-p and 
n-n momentum correlation functions (solid symbols) and the 
anti-one (open symbols) are shown in each panel. 


Figure 3 shows the centrality and system-size depen- 
dences of the momentum correlation functions for light 
(anti)cluster in similar condition as in Fig. 2. Figure 3 
(a) and (c) present the momentum correlation functions 


of d-d (d-d) and t-t for 19" Au +237 Au collisions at dif- 
ferent centralities of 0 — 10 96, 10 — 20 96, 20 — 40 96, 
40 — 60 %, and 60 — 80 % at \/s~n = 39 GeV. The d-d 
(d-d) momentum correlation functions exhibit less than 
unity at lower relative momentum q in Fig. 3 (a) and (b), 
which is caused by the Coulomb repulsion. Then, with 
increasing relative momentum q, the anti-correlation be- 
tween two-(anti)deuteron pairs is more complex. The 
two-triton momentum correlation functions are less than 
unity with increasing relative momentum q as shown in 
Fig. 2 (c) and (d), which is caused by only the Coulomb 
potential in the Lednický and Lyuboshitz code [45-47]. 
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FIG. 3. Similar to Fig. 2 but for the light (anti)cluster pairs. 
(a) (b) d— d momentum correlation functions (solid symbols) 


and the anti-one (open symbols), (c) (d) t — t momentum 
correlation functions (solid symbols). 
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FIG. 4. Similar to Fig. 2 but for the nonidentical particle 


pairs. (a) (b) n-p momentum correlation functions (solid sym- 
bols) and the anti-one (open symbols), (c) (d) p-p momentum 
correlation functions (solid symbols). 


The antideuteron—antideuteron momentum correlation 
function also shows an exact similar shape with deuteron 
pairs due to the similar phase-space distributions be- 
tween deuteron and antideuteron. Due to significant 
less yields of tritons which induce too large error, the 
antitriton—antitriton momentum correlation function is 
not shown in the present work, which should be observed 
as the same trend with triton pairs. Fig. 3 (a) and (c) also 
compare five centralities of 0— 10 96, 10—20 96, 20— 40 96, 
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FIG. 5. 


Similar to Fig. 4 but for nonidentical light (anti)nuclei. (a) (b) p-d momentum correlation functions (solid symbols) 


and the anti-one (open symbols). (c) (d) p-t momentum correlation functions (solid symbols), (e) (f) d-t momentum correlation 


functions (solid symbols). 


40 — 60 96, and 60 — 80 96 for the momentum correlation 
functions of two light (anti)clusters. The more suppres- 
sion of the d-d (d-d) and t-t correlation functions is clearly 
visible in peripheral selection. These results also indicate 
that light (anti)cluster emission occurs from a source of 
smaller space—time extent for peripheral collision, which 
is the similar to Fig. 2 (a) and (c). In Fig. 3 (b) and (d), 
an enhancement strength of the momentum correlation 
function for d-d (d-d) and t-t is also observed when light 
(anti)cluster pairs emitted from smaller systems. How- 
ever, for small systems such as B and O, the sensitivity 
seems disappear. 


C. Nonidentical light (anti)nuclei momentum 
correlation functions gated on centrality and 
system-size 


Now we investigate centrality and system-size depen- 
dence of the nonidentical (anti) particle momentum corre- 
lation functions, such as n-p (n-p), p-p, p-d (p-d), p-t and 
d-t. Fig. 4 (a) and (c) show results for the momentum 
correlation functions of n-p (n-p) and p-p for the same 
centrality classes as Fig. 2. The same centrality depen- 
dence is also clearly seen in Fig. 4 (a) and (c). Because 
of the strong attractive final state interaction between 
n and p, the n-p (n-p) momentum correlation functions 
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FIG. 6. The velocity-gated momentum correlation func- 
tions in relative momentum (left) and velocity difference (Av) 
spectra (right) for n-p and p-p as a function of five different 
centralities at mid-rapidity (|y| < 0.5) in /sww = 39 GeV 
79. Au --19/ Au collision. The velocity conditions are indicated 
in each panel, Av > 0 (solid symbols) and the anti-one (open 


symbols). 
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FIG. 7. The velocity-gated momentum correlation functions 
in relative momentum (left) and velocity difference (Av) spec- 
tra (right) of n-p and p-p pairs for 0 — 10 % central collisions 
of 19B -I? B, 190 4-19 O, $8Ca +59 Ca as well as 12" Au --12" Au 
systems at J/sww = 39 GeV. The velocity conditions are indi- 
cated in each panel, Av > 0 (solid symbols) and the anti-one 
(open symbols). 


show a strong positive correlation for small values of the 
relative momentum q in Fig. 4 (a) and (b). Fig. 4 (c) 
shows results for proton—antiproton momentum correla- 
tion functions, which are different from the results for 
proton pairs in Fig. 2 (c), however, qualitatively agrees 
to the experimental results in Ref. |78, 79]. In addition, 
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FIG. 8. Ratios of the velocity-gated momentum correlation 
functions in relative momentum (left) of n-p (a) and p-p (c) 
pairs as a function of five different centralities at mid-rapidity 
(|y| < 0.5) in /sww = 39 GeV 19" Au 4-H$' Au collision. Ra- 
tios of the velocity-gated momentum correlation functions in 
relative momentum (right) of n-p (b) and p-p (d) pairs for 0 
— 10 96 central collisions of 19B --1? B, 190 4-19 O, 58Ca 4-20 Ca 
as well as 19 Au +72" Au systems at V 5NN = 39 GeV. 


Fig. 4 (b) and (d) show system-size dependence of n- 
p (n-p) and p-p momentum correlation functions, which 
is almost unanimously with the identical (anti)particle 
one in Fig. 2 (b) and (d). We can also observe an en- 
hancement strength of particle pairs momentum corre- 
lation function in smaller systems. In the same way, 
we also investigate the effect of different centralities and 
system-size on the momentum correlation functions of 
nonidentical light (anti)nuclei. The p-d (p-d), p-t and 
d-t momentum correlation functions in Fig. 5 (a), (c) 
and (e) are all characterized by an anti-correlation. For 
the p-d (p-d) momentum correlation functions in Fig. 5 
(a), the anti-correlation shape is a little unlike to the 
proton—deuteron momentum correlation function in the 
intermediate energy heavy-ion collision [34, 35], indicat- 
ing that a competition between the s-wave attraction and 
the Coulomb repulsion. The correlation functions of p-t 
and d-t in Fig. 4 (c) and (e) also display the trend of be- 
low unity due to the dominant Coulomb repulsion, which 
is similar to the previous results in intermediate energy 
heavy-ion collisions |34, 35]. In Fig. 5 (b), the system-size 
dependence of p-d (p-d) momentum correlation functions 
is shown, while we can observe an enhancement of p-d 
(p-d) momentum correlation function in smaller systems. 
In Fig. 5 (d) and (f), the p-t and d-t momentum correla- 
tion functions appear more sensitive to system-size only 
in the large system such as Au and Ca. 
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FIG. 9. Similar to Fig. 6 but for p-d (a) (b), p-t (c) (d) and d-t (e) (f) pairs. 


D. Velocity selected nonidentical light nuclei 
momentum correlation functions 


The momentum correlation functions of unlike par- 
ticles can provide an independent constrain on their 
mean emission order by simply making velocity selec- 
tions [22, 32, 33, 80, 81]. The principle of comparing the 
velocity-gated momentum correlation functions for the 
nonidentical particle pair to infer their average emission 
order is as follows. Here the two nonidentical particles 
are named by “a” and “b” , respectively. If the ve- 
locity of “a” particle is lower than “b” particle, the 
(anti)correlation will be stronger when the "a" particle 
is emitted averagely early than the "b" particle, be- 
cause the space-size between them is reduced during the 
flight and the final-state interaction (FSI) is enhanced, 
and vice versa. In addition, the velocity difference (Av) 
spectrum between the two nonidentical particles is also 
sensitive to the mean emission order. Fig. 6 presents the 


velocity-gated momentum correlation functions and ve- 
locity difference (Av) spectra of unlike particles pairs n-p 
and p-p for 19" Au +23” Au collisions at different central- 
ities of 0 — 10 76, 10 — 20 96, 20 — 40 96, 40 — 60 %, and 
60 — 80 96 at \/sww = 39 GeV. In Fig. 6 (a) and (c), the 
centrality dependence on the velocity-gated momentum 
correlation functions of n-p and p-p are similar to Fig. 4. 
In Fig. 6 (a), the momentum correlation function for n-p 
pair with v, > vp is similar to one with the reverse situ- 
ation. The symmetry of velocity difference (Av) spectra 
for n-p pairs is shown in Fig. 6 (b). The results demon- 
strates that the average emission sequence of neutrons 
and protons are almost the same and have no sensitive 
to the centrality. In Fig. 6 (c), the momentum correla- 
tion function for p-p pair with vp > vy is slight higher 
than one with the reverse situation. ‘The slight asym- 
metry of velocity difference (Av) spectra for p-p pairs is 
shown in Fig. 6 (d). The corresponding relation indicates 
that the mean order of emission sequence between pro- 
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FIG. 10. Similar to Fig. 7 but for p-d (a) (b), p-t (c) (d) and d-t (e) (f) pairs. 


ton and antiproton may be a little different and have also 
no sensitive to the centrality. In Fig. 7 (a), the momen- 
tum correlation functions for n-p pairs with v,, > vp are 
always similar to one with the reverse situation with in- 
creasing system-size. The symmetry of velocity difference 
(Av) spectra for n-p pairs in different systems is shown in 
Fig. 7 (b). The comparison of velocity-gated momentum 
correlation functions indicates that the average emission 
sequence between neutrons and protons can be always 
identical for different centrality and system-size, which 
are also shown by their ratios in Fig. 8 (a) and (b). In 
Fig. 7 (c) and (d), the comparison of velocity-gated mo- 
mentum correlation functions for p-p indicates that the 
mean order of emission sequence between protons and an- 
tiprotons may be a little different and has no dependence 
on system-size, which are also shown by their ratios in 
Fig. 8 (c) and (d). 


Fig. 9 and Fig. 10 also show centrality and system- 


size dependence of velocity-gated momentum correlation 
functions and velocity difference (Av) spectra of p-d, p- 
t and d-t pairs, respectively. For p-d and p-t pairs, the 
momentum correlation functions with uy < va (Up < vi) 
is stronger than the one with the reverse situation vp 
> va (Up > ve) in Fig. 9. The comparison of the two 
anti-correlation strengths gives that the mean order of 
emission that protons are emitted averagely earlier than 
deuterons and tritons according to the above criteria. 
The similar trend for d-t pairs is not obviously, where the 
momentum correlation function with vg < v, is stronger 
and deuterons are emitted averagely earlier than tritons 
especially in peripheral collision. However, the average 
emission sequence of protons, deuterons, and tritons is 
opposite to the emission order in previous results for the 
intermediate energy heavy-ion collisions |11, 34, 35, 80]. 
At the same time, Fig. 9 presents velocity difference spec- 
tra for p-d, p-t and d-t pairs, respectively. The velocity 
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FIG. 11. 


difference spectra are all asymmetric due to the mean 
emission order. In addition, an enhancement difference 
between the momentum correlation functions for p-d (p-t 
or d-t) pairs with vp > vq (Up > v, or uq > v) and ones on 
the reverse situation with larger centrality, which man- 
ifests the larger interval of the mean emission order for 
unlike light nuclei in peripheral collisions. Their ratios in 
Fig. 11 (a), (c) and (e) can also illustrate the above phe- 
nomenon. The system-size dependence for p-d, p-t and 
d-t pairs can be found that the momentum correlation 
function with vy < vq (vy < v or Va < Uz) is stronger 
than the one with the reverse situation vp > vq (Up > t 
or vq > v+) in Fig. 10. Correspondingly, the velocity 
difference spectra for p-d, p-t and d-t pairs are all asym- 
metric about Av = 0 caused by the average emission or- 
der in Fig. 10. Therefore, protons are emitted averagely 
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(b) 39GeV 0-10% Av =v V, 
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Similar to Fig. 8 but for p-d (a) (b), p-t (c) (d) and d-t (e) (f) pairs. 


earliest and deuterons are emitted averagely earlier than 
tritons in smaller system-size collision. ‘The system-size 
dependence of the velocity-gated momentum correlation 
functions is also clearly seen by their ratios in Fig. 11. 
With decreasing system-size, we can also observe an en- 
hancement difference between the momentum correlation 
functions for p-d (p-t or d-t) pair with vp > va (Up > vi 
Or Ug > 14) and one with the reverse situation in Fig. 11 


(b), (d) and (f). 


IV. SUMMARY 


In summary, with the AMPT model complemented 
by the Lednicky and Lyuboshitz analytical method, we 
have constructed the momentum correlation functions of 


light (anti)nuclei formed by the coalescence mechanism 
of (anti)nucleons at different centrality and system-size 
in snn = 39 GeV heavy-ion collisions. We present 
a comparison of proton— proton and proton—antiproton 
momentum correlation functions with the experimental 
data from the RHIC-STAR collaboration [78, 79]. Tak- 
ing the same transverse momentum and rapidity phase 
space coverage corresponding to the experimental situa- 
tion as well as the maximum hadronic rescattering time 
of 700 fm/c in AMPT, it is found that the p-p and 
p-p momentum correlation functions simulated by the 
present model can match the experimental data. We 
further study centrality and system-size dependence of 
momentum correlation functions for identical and non- 
identical light (anti)nuclei pairs, respectively, which is 
in the condition of the maximum hadronic rescattering 
time of 100 fm/c in AMPT. The shape of momentum 
correlation functions for light (anti)nuclei pairs is consis- 
tent with previous works [13, 28, 34, 35, 78, 79], which 
is caused by both QS and FSI. The similar structure be- 
tween light nuclei momentum correlation functions and 
anti-ones indicates that the interaction between them 
are the same, which has been confirmed in Ref. [28] 
only about proton and antiproton. The centrality de- 
pendence of momentum correlation functions for light 
(anti)nuclei is investigated by #2’Au +737 Au collisions 
at different five centralities of 0 — 10 96, 10 — 20 %, 
20 — 40 96, 40 — 60 96, and 60 — 80 % at \/swn = 39 GeV. 
It is found that with increasing centralities from center 
to peripheral, the momentum correlation functions for 
light (anti)nuclei become stronger, which are probably 
emitted from smaller source. The momentum correla- 
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